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Abstract: The next generation of ultra high energy cosmic ray experiments will 
probably detect several dozens of events clustered around the direction towards each 
of the most powerful extragalactic sources. We develop a method which could make 
possible to reconstruct, from the arrival directions and energies of the clustered 
events, the strength and coherence properties of the magnetic field along the line 
of sight towards the sources. The method exploits peculiar signatures arising from 
magnetic lensing effects, such as the strong flux magnification of multiple images 
around caustics. We also discuss how to obtain information about the cosmic ray 
composition, and apply this method to samples of simulated data. 
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1. Introduction 

The origin and nature of ultra high energy cosmic rays (UHECRs), after several 
decades of study, is still puzzling |l|, H . If they are charged particles, such as protons 
or heavier nuclei, their trajectories are affected by the magnetic fields along their 
path from the sources to the observer. For energies below Z x 10^^ eV (Ze is the CR 
electric charge) they are confined by Galactic magnetic fields. For larger energies 
the galactic field is not strong enough to confine them, which combined with the 
lack of any significant excess from the galactic disk suggests that UHECRs are most 
probably of extragalactic origin. However, due to the effects of the magnetic fields 
present along their path, they do not point to their birth places, complicating the 
task of identifying their sources. The small angular scale clustering observed in the 
AGAS A data already hints to the existence of UHECR point sources ||, ^ , but 
the small present statistics does not allow to solve the source identification question. 
With the next generation of detectors, like AUGER and EUSO larger clusters 
are expected, having each several dozens of events. If this were the case, it would not 
only become possible to reconstruct the source positions, but also the data could be 
used to obtain information about the magnetic field along the CR path |10|, |TT], [12 



In this paper we discuss some strategies which can be adopted to reconstruct the 
main parameters of Galactic magnetic fields from a set of events originating from 
an extragalactic point source, and test them using simulated data. The knowledge 
of the magnetic field responsible for the deflections allows to better reconstruct the 
locations of the sources and, moreover, can also be helpful to do spectroscopy of CRs, 
making possible the measurement of their charges. The method we develop here 
profits from several specific features imprinted on the events by the magnetic fields, 
which not only deflect the trajectories but also lead to strong lensing phenomena, 
including the formation of multiple images and energy dependent magnifications or 
demagnifications of the CR fluxes |TB|, |T3| . 
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2. Magnetic lensing 



The magnetic field in the Galaxy is known to have a regular and a turbulent com- 



ponent [jT6l. The regular field in the disk follows the spiral arms structure, with 
reversals taking place from arm to arm. The local value is Breg — 2 nG. The possible 
existence of an extended halo field (of few kpc scale height) is not settled down yet. 
The typical defiection in the direction of propagation of a CR particle of charge Ze 
and energy E produced by a regular magnetic field is 

' ds B 



40EeV 



3 kpc 2 yuG 



(2.1) 



where 1 EeV = 10^^ eV. 

The turbulent galactic magnetic field component has a root mean square am- 
plitude larger than the regular field, Brms — 1^3i?r-eg, and the largest turbulence 
scale is Lmax — 100 pc, while the minimum one is believed to be much smaller 
{Lmin <^ Lmax)- It IS ofteu modeled as a Gaussian random field with a Kolmogorov 
spectrum |17] (i.e. such that the energy density satisfies dE/dk oc k~^^^). The mean 



defiection of particles moving in the turbulent field vanishes, while the root mean 
square defiection is given by [|1^ 



5_.1.4°^^J^W^, (2.2) 
E/Z 4/iGy Skpcy 50pc' ^ ' 

where Lc is the coherence length. For a narrow band {Lmin — Lmax) spectrum, one 
has Lc — Lmax/'^i while for a broad band {Lmax ^ Lmin) Kolmogorov spectrum it is 

Lc — Lmax / ^ ■ 

We have written the random and regular field defiections in terms of the charac- 
teristic values of the Galactic fields, thus the numerical value of the defiections are 
typical expected values. We see that the regular field generally produces the domi- 
nant effect on the deflection, which leads to a coherent displacement of the apparent 
position of the source that is inversely proportional to the energy. Compared to this, 
the random fleld produces rather small changes in the arrival direction of the CRs. 
However, as the charged particles propagate through the magnetic fleld, they are 
not only deflected, but their flux is also focused or defocused due to the differen- 



tial deflections of neighboring paths O, 14 1, with the effects being larger at smaller 



energies. As lower energies are considered, CR particles from one source can have 
more than one path leading to the detector, each traversing an uncorrelated magnetic 
field patch, and this leads to the appearance of multiple images of the source. This 
typically occurs when 6rms — Lc/L |T3|, what allows to define a critical energy Ec, 
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around and below which the formation of multiple images is very likely, through^ 



hj L 



Its numerical value is given by^ 



E^^Zm EeV^ ( ,1^. (2.4) 



4 /iG V3 kpcy V L 



-'c 



The new images appear in pairs, at an angular distance ~ L^jL from the position 
of the original source image and with a large magnification of their fluxes. The peak 
in the amplification when new images appear leads to small angular scale clustering 



and to an excess of events in an energy bin close to Ec [15 . 

The overall picture can be summarized as follows: the magnetic fields present 
along the CR trajectories change the mapping between the observed arrival directions 
and the source ones. At the highest energies, the effects are small and hence the 
mapping is close to the identity. Going down in energy, deflections arise due to 
the turbulent and the regular magnetic field components. The turbulent component 
leads to small rms deflections, but it gives rise to a network of lines in the sky around 
which the fluxes of potential sources would be significantly magnified. This network 
corresponds to the directions where folds in the mapping (caustics) develop. The 
caustics appear at energies near (and below) 2 £'c, giving rise to multiple images of the 
sources lying close to those directions . The regular field, which is responsible for 
larger deflections, leads to a global displacement of this network of caustics relative 
to the source directions. As a result of this drift caused by the regular field, a 
given CR source will experience successive magnifications and demagnifications of 
its flux at different energies caused by the motion of the caustic network across the 
sky. For decreasing energies, the continuous appearance of new images also leads to 
additional peaks in the spectrum, and strong lensing effects associated to the regular 
field component will also show-up. 



3. Magnetic field reconstruction 



The properties of the magnetic deflections and lensing effects can be used to de- 
termine the magnetic field parameters from the study of clustered UHECR events 
having their origin in isolated point-like sources. For this purpose one should start 



^The regular field can also lead to multiple images, with its characteristic scale of homogeneity 
(few kpc) playing the role of the coherence length. 

^If the amplitude of the turbulent field were not constant, but had instead some modulation 

^\fTj in eq. ( |2.2| ) should be replaced by | /q^ ds i3,^„^(s)|^/^, and 



along the Galaxy, the factor B, 
BrmsLi^^^ in eq. (2^) should be replaced by 3 Ax x ds B^^g{s) 



1/2 
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by defining the regions containing the excess events (which have to be larger than 
at least a few times the angular resolution of the detector but yet not too large to 
avoid excessive background events coming from other independent nearby sources). 
Different strategies must be considered, depending on the relative strength of the ef- 
fects due to regular and random components, and in what follows we develop several 
ideas to perform this reconstruction. 

Consider first the hkely case in which the defiection caused by a regular com- 
ponent (either galactic or extragalactic) is larger than the rms defiection imprinted 
by the turbulent galactic field. In this case, the overall angular motion of the source 
images as a function of energy can be used to reconstruct the integral along the line 
of sight to the source of the perpendicular component of the regular magnetic field. 
In order to do this, one first finds the overall direction of motion of the images caused 
by the regular field by fitting a straight line {ay — a + bax) to the event coordinates 
ax,y The integral along the line of sight of the perpendicular component of the 
regular magnetic field will just be at a right angle with respect to this line, and its 
absolute value can be obtained choosing now the coordinates ay (in the direction 
along the deflections produced by the regular field) and a± (orthogonal to it), and 
plotting a\\ vs. 1/E for the set of selected events. One should have, ignoring at this 
step the random field component, that 

^. 

ail . ~ ao + TT^-L (3.1) 
' Ei 

with 

ds B'^^s 



X, = ±|e 



/" ds X B^^^l ~ ±1.3°(40 EeV) [ 
Jo Jo 



X 



kpc /.G • ^^-^^ 

If the CRs have all identical electric charge Z, the events would fit nicely to 
just one straight line, a\\ = + c/E. The slope c would give the value of Z K±, 
which essentially provides the magnitude and direction of the regular magnetic field 
orthogonal to the l.o.s. multiplied by its relevant scale length, modulo the CR electric 
charge. The value of ao just gives the original location of the source {a\\ — ckq, 
a± — 0). After the source location is estimated, it proves convenient to change 
variables to sign{K±) (ay — ao), defined such that it tends to zero at large 

energies and increases as E decreases. In what follows we will use the name ay just 
to refer to this shifted variable. 

We illustrate this procedure using simulated data, obtained through a ray shoot- 
ing technique, in which a large number of antiparticles were thrown isotropically 
from Earth, their trajectories numerically integrated across a distance L within a 
homogeneous magnetic field plus a turbulent component, and were finally recorded 
only if they pointed within a small cone (of aperture ~ 5rms/10) from the direction 
to a fixed 'source'. The source was assumed to inject CRs with a differential energy 
spectrum dN/dE oc E~'^-'^ . About 100 simulated events with energy above 40 EeV 
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were selected, assuming the detector efficiency to be energy-independent. We dis- 
play the 'theoretical' results from the simulation (with black dots) as well as more 
realistic ones (open circles) obtained by adding noise to the energy (corresponding to 
a detector energy resolution of 10%) and to the angles (corresponding to an angular 
resolution of 0.2° in each direction). 

In the case illustrated in Figure |I| the homogeneous magnetic field had strength 
2 yuG along the direction, the turbulent component had -Brms = 4 with just 
one turbulence scale L^ax = Lmm = 100 pc (and thus Lc = 50 pc), and the distance 
traversed by the CRs (assumed to be protons) within these fields was L = 3 kpc. The 
source actual position was in the direction (a^. = 0, = 0). The first panel displays 
the simulated data, with the straight line being the reconstructed overall direction 
of motion of the images and the star being the reconstructed source location. The 
second panel plots instead a\\ vs. 1/E, while the third one will be discussed farther 
below. 




Figure 1: First Panel: Angular distribution of the simulated events for a regular field 
Breg = 2 /xG, a turbulent field with Brms = 4 /xG, L^ax = -^^mm = 100 pc and L = 3 kpc. 
White circles have a Gaussian noise with dispersion 0.2° added in each direction, as well 
as a 10 % uncertainty in the energy, to simulate the detector sensitivity, while black circles 
have no noise added. The size of the circle grows with the energy of the event. The asterisk 
indicates the reconstructed position of the source, and the solid line the fit to the deflection 
due to the regular field. Second panel: ay vs. 1/E [EeV~^] and the linear fit to the data. 
Third panel: Number of events vs. aj| divided in ten angular bins. 

In Figure |^ we display the results for another simulation with the same parame- 
ters for the regular field and same root mean square amplitude of the turbulent field 
but with a Kolmogorov spectrum with Lmax = ^OLmm = 100 pc (corresponding to 
Lc = 25 pc). The result for the best fit of the numerical reconstruction was for both 



- 5 - 



cases Z| ds X B| ^ 6.1 fiG kpc, in good agreement with the values used in the 
simulations. 

In the case of a mixed composition, the fit of a\\ vs. 1/E to just one straight line 
should fail (excessive x^), since one expects to see the superposition of several straight 
lines, all with the same value of but different slopes, with ratios fixed by the ratio 
of the different electric charges Zj involved^. In this case one can separate the events 
in groups, above and beyond limiting slopes, and then fit each group independently 
to an expression of the form given in Eq. (|3.1| ). This allows to obtain information 
on both K± and the ratio of charges Zi involved, besides the relative amounts of 
different nuclei from the relative number of events on the different lines. Notice that 
the reconstructed fraction of different nuclei reflects the composition upon arrival to 
Earth, which may differ from the fraction injected by the source within the same 
energy range because the flux magnification due to magnetic lensing depends on the 
combination E/Z. Thus, the flux of different nuclei is magnified by a different factor 
at a given energy [H. 
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Figure 2: Same as Figure || for a turbulent field with a Kolmog orov spectrum and Lmax — 

WLmin = 100 pC. 

In Figure ^ we exemplify this procedure for the same values of the fields and 
Kolmogorov spectrum with 60% of protons and 40% of Lithium {Z = 3). (The actual 
fraction of Lithium injected by the source was 50% in this case.) The two slopes in 
a\\ vs. 1/E are clearly separated and the reconstructed values are Zi\ ds x B| ^ 
5.9 yuG kpc, and Z2I Jq ds x B| ^ 18.6 jjG kpc, implying that Z2/Z1 3.1. 

In Figure ^ we show the analogous results for 80% of protons and 20% of Helium, 
for which we obtain Zi\ ds x B | ^ 6.2 /xG kpc, and Z2I /q^ ds x B| ^ 12.3 /xG kpc, 
^The value of in the case of mixed composition was just determined by eye, looking for the 
value of a|| towards which the different lines converged. 



- 6 - 



implying that ^2/^1 ~ 2. The reconstructed direction of the regular magnetic field 
was off by only a few degrees in all cases. Let us mention that we performed several 
other simulations with different realizations of the random field component, obtaining 
reconstructed values of similar quality. 




a, 1/E a. 



Figure 4: Same as Figure ^ with mixed composition, 80% of protons and 20% of Helium. 

Having determined the overall motion of the images due to the regular magnetic 
field component, we turn to analyze the effects of the turbulent field. At high energies 
the main effect is the formation of multiple images, with their associated peaks 
manifesting in the observed flux. The critical energy of this process, E^. (given by 
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eq. (|2.4| )), can be determined from the location of the first lensing peak, which is 
expected to he at energies somewhere between 2Ec and Ec/2. We then take bins in 
a\\ covering the range of angles where clustered events are observed above a certain 
energy threshold. The bins should be larger than the angular resolution of the 
detector and such that a sizeable number of events per bin result (taking e.g. ten 
angular bins if the number of clustered events is of the order of one hundred). Now 
one can plot the number of events in each bin vs. a\\. If there were only deflections 
but no flux (de)magnification due to lensing effects, the regular field would just lead 
to AiV/Aa|| oc ap where P characterizes the slope of the CR energy spectrum, with 
dN/dE oc 

The advantage of plotting AA^ in terms of a\\ (and not e.g. vs. E"^) is that 
regardless of the CR composition the first lensing peak always appears at the same 
value of a\\. Referring to this value as ajj"^"'^, one has then (up to an overall factor of 
two) 

Er. \K<\ 



^ ^peak 



(3.3) 



The third panels of Figures illustrate the results of this method applied to 
the simulated data. In Figure |l| the first lensing peak associated to multiple image 
formation appears at ay*^"*^ ~ 6.4° and corresponds to the sharp peak. It leads to a 
critical energy value E^/ Z ^ 49 EeV, in good agreement with the parameters used in 
the simulation (which correspond to Ec ~ 60 EeV). A shallower and more symmetric 
excess of events is also evident at smaller deflection angles, a^^^"'^ ~ 2.5°, and it is due 
to a high magnification region crossing the source location at high energies (larger 
than ^ 2Ec, when folds were not yet formed). 

In Figure |^ the first peak appears at aji'^"'^ ~ 4.3°, corresponding to Ec/Z 
72 EeV (while the simulation parameters lead to Ec/Z ^ 85 EeV). For the mixed 
composition cases (Figures |^ and |) the corresponding values are Ec/Z ^ 85 EeV 
and Ec/Z ^ 87 EeV respectively. 

The angular distribution of the events can also be used to estimate Lc/ L [L here 
is the distance traversed along the turbulent component, which needs not be the 
same as that across the regular one). There are two independent ways to do this. 
One is based on the fact that at energies ~ Ec the typical angular separation in the 
network of caustics of the turbulent magnetic field is 4Lc/L JT^. Since the effect of 
the regular field induced defiections is to move these caustics, relative to the source 
direction, the typical separation in a\\ between consecutive lensing peaks is also equal 
to ALc/L. If at least two lensing peaks are observed in the data from one cluster, 
then Lc/ L can be estimated. This is the case in the simulations presented above. In 
the case with just one wavelength depicted in Figure |] (which has Lc = 50 pc), the 
separation between peaks is of order 4°, as it can be seen in the third panel, leading 
to Lc/L ^ 0.017. In the Kolmogorov spectrum cases (for which Lc — 25 pc) depicted 
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in Figures ^, ^ and H, the separation among peaks is between 2° and 3°, leading to 
Lc/L values in the range 0.008-0.013. All these values are within 50% of the actual 
values used in the simulations. 

Notice that it is likely to observe the first lensing peak rather close to 2i?c5 and to 
observe more than one lensing peak in the range of energies between 2Ec and -E'c/2, 
if the deflection caused by the regular magnetic field component is strong enough. 
Indeed, the network of caustics is already well defined at energies around 2Ec. Al- 
though only a small percentage of all potential source directions are sufficiently close 
to a caustic at energies around 2Ec to have their flux noticeably magnified, if the 
deflection caused by the regular magnetic field component is significant, a caustic 
will reach the source position at energies not far below 2Ei.. For instance, if the reg- 
ular field deflection is at least 3 times larger than the rms deflection caused by the 
random field, then the source apparent position is displaced, as the energy decreases 
from 2Ef. down to -Ec/2, by more than ALc/L. Since this is the typical separation 
between caustics, it is likely to produce more than one lensing peak in this energy 
range. Notice however that if the deflection by the regular field were too strong, the 
caustics may be crossed too fast, leading to a small integrated effect and thus a less 
noticeable peak. 

Alternatively, L^/L can be determined from the fact that it is just the typical 
angular scale characterizing the separation between the images present in the first 
lensing peak. In this peak there are in principle one image corresponding to the 
principal image which is also present at higher energies, and two new images which 
are both in the same location of the sky and have appeared separated from the 
previous one by an angular scale of order Lc/L. Since this angular separation may 
be quite small, this approach can be used only when the angular uncertainty in the 
events is smaller than this scale. For the example shown in Figure |^ this separation is 
about 1°, leading to an Lc/L estimate similar to the one obtained with the previous 
method. For the Kolmogorov spectrum cases, that leads to smaller Lc/L, we see 
from the first panel in Figures 0, ^ and |^ that to apply this method in cases like 
these would require a quite precise angular resolution, probably better than the one 
achievable with AUGER (~ 0.3°), but possibly within the aim of some proposed 
detectors |T^. 

Once Ec/Z and Lc/L have been estimated, we can combine them to obtain an 
estimate of Brms L using Eq. (^.41). In our examples, the recovered values fall in the 
range Brms L ~ 10-13 fiG kpc, in good agreement with the value 12 fiG kpc used in 
the simulations. 

If the deflection caused by the regular magnetic field is smaller than the rms 
deflection caused by the turbulent field at the same energy, the overall motion of 
the images does not provide reliable information about the regular magnetic field 
component, other than an upper bound on its integrated effect. However, it might 
still be possible to better reconstruct the source location by looking at the highest 
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energy events, i.e. at energies above those at which the multiple images appear. 
Indeed, as long as the deflections are small {6 -C Lc/L), the random field component 
still produces a coherent defiection of the images proportional to {E/Z)~^, and hence 
unless there is a significant mixture of different CR compositions there will be a clear 
signal allowing to reconstruct the source location. One could also extract information 
about the turbulent field through the spectral features imprinted by the lensing peaks, 
and by the angular clustering properties of the events. The critical energy can be 
estimated as the highest energy at which a peak is observed in the spectrum. With 
sufficient angular resolution, Lc/L can be estimated from the angular separation 
between the images present in the first lensing peak, while information about the 
composition can be extracted from the observed energy ratios of the events clustered 
in the peaks ]15 . 




Figure 5: Simulated data in a regular field Breg = 0.1 ^G, and other parameters as in 
Fig. |l|. The injected composition is 100% protons (left panel), and a mixture with 40% 
Beryllium (right panel). 

These general features are illustrated in Figure ^ which corresponds to a ray 
shooting simulation with a small orthogonal component of the regular magnetic field 
(leading to J B^^^ds = 0.3 fiG kpc). The random component has Bms = 4 /iG (and 
spectrum with L^m = L^ax = 100 pc) and L = 3 kpc. Events with energies larger 
than 40 EeV and injection spectrum oc are displayed. The left panel is for 

CR protons alone. The highest energy events move away from the source location 
(here the origin) with a defiection oc E~^, and at energies near 60 EeV two pairs of 
secondary images appear. They are displaced by ~ 1.5° from the principal image, 
hence leading to an estimated L^/L ~ 0.026, in good agreement with the input value 
0.017. For the right panel, which has a mixed composition of protons and 40% Be at 



- 10 - 



injection, other secondary images appear for ii^ > 40 EeV, dispersed over a region of 
~ 3°, and they contain exclusively the heavier Be nuclei, since the protons in them 
have energies below 40 EeV. In general each of the secondary images has a narrow 
rigidity distribution due to the lensing origin of the peaks (and hence if two different 
elements were contributing to an image in the considered energy range, the events 
would split into groups with energy ratios given by the ratios of their charges). 



4. Conclusions 

One of the most exciting perspectives for the Auger Observatory currently under 
construction, as well as for other future experiments, is that they may inaugurate 
the era of UHECR astronomy. A hint that this may be the case is the small angu- 
lar scale clustering observed by AG AS A, which may be the first indication that a 
few UHECR sources are significantly more powerful than the otherwise apparently 
isotropic background. 

UHECR astronomy is inevitably tied to the magnetic fields along the line of sight 
towards the most powerful sources, if the bulk of their emission is in the form of 
charged particles, such as protons or heavier nuclei. Conservative estimates of galac- 
tic and intergalactic magnetic fields indicate that the deflection of protons should 
be small at the highest energies, around and above 10^° eV, so that they should 
truly point to their birthplace. However, even protons with energies around 10^° eV 
may have their flux significantly (de)magnified by lensing effects, for instance by the 
turbulent component of the magnetic field in the Milky Way. A heavier component 
of UHECRs would undergo deflections larger than protons at comparable energies, 
and would be subject to similar flux (de)magnification effects at higher energies. 

Consequently, the behaviour of the transition towards UHECR astronomy at the 
highest energies is a very rich source of information both on the source properties as 
well as on the characteristics of the intervening magnetic fields. This has been ex- 
ploited, for instance in |ll| , as a tool for the reconstruction of intergalactic magnetic 



field parameters, particularly in the case of bursting sources of UHECRs , where 
the energy-dependent time delays provide additional handles to perform the task. 

Here we developed a strategy to reconstruct the parameters of intervening mag- 
netic fields that profits from the strong lensing effects that take place at energies 
of the order of Ec, given by Eq. ( p^.4| ), around which it is likely to observe strongly 
magnified multiple images of single UHECR sources. We have shown that it is likely 
that a sort of 'magnetic spectroscopy' of UHECR sources may allow to measure the 
strength of the magnetic field along the line of sight, its coherence properties, and 
the relative amounts of nuclei with different electric charges among the observed 
CRs. Our analysis assumed the sources to be steady, what means here essentially 
that the typical time delays between different images are smaller than the timescales 
of emission of the sources. 
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The method works if the number, strength and location of sources, the UHECR 
composition, and the intervening magnetic field parameters are such that clusters 
with a large number (several dozens) of events can be identified, without excessive 
background from events with origin in nearby independent sources. The angular and 
energy resolution required is also dependent on the actual values of the parameters 
involved. We have shown, through illustrative examples based in simulated data, 
that the method is likely to be applicable for realistic values of the parameters of 
the galactic magnetic field. The overall angular motion of the events as a function of 
energy allows to reconstruct the integral along the line of sight of the perpendicular 
component of the regular magnetic field, modulo the CR electric charge. If this 
can be done for several sources, one may then even map the strength and extent 
of the magnetic field along different directions, and this could help to establish its 
overall distribution.^ The reconstructed magnetic field also provides a handle to 
measure the relative number of CRs with different electric charges. A 'spectroscopic' 
analysis of the peaks in the number of events as a function of this overall motion 
(or as a function of energy when the effect of the regular component is not strong 
enough) provides a measure of the critical energy E(./Z at which strong lensing 
phenomena occur. The enhanced angular clustering of the events at the energies 
where strong flux magnification takes place can also be used to estimate Lc/L, the 
ratio between the coherence length of the turbulent magnetic field component and 
the path length traversed by the CRs across it. Ec/Z and Lc/L provide a measure 
of BrmsL, independent of the CR electric charge. The fact that the magnification 
peaks due to magnetic lensing occur at energies which are in direct proportion to the 
electric charge Z provides a second handle to grasp the CR composition. It induces 
a strong correlation of arrival directions at energy ratios fixed by the ratios of the 
CR electric charges. 

Potentially, this method could also extract information about the different scales 
involved in the turbulence of the intervening magnetic field, if the experiments had 
sufficient energy and angular resolution (and, of course, statistics). If the spectrum 
is broad band, long and short wavelengths lead to amplification peaks at somewhat 
different values of E/Z |T^. The density of lensing peaks could thus be used to 
obtain information on the magnetic field spectral properties. However, in realistic 
experiments the very narrow lensing peaks due to short wavelengths are likely to go 
unnoticed. 

The method may also provide information on extragalactic magnetic fields. No- 
tice that for an extragalactic turbulent component one has an associated critical 



^Notice that the information reconstructed with UHECRs is complementary to that provided 
by Faraday rotation measurements, which are sensitive to the magnetic field component parallel to 
the line of sight. 
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energy 

19 ^ Brms f L /MpC 



E..Z2xlO'»EeV^^_J (4.1) 

Hence, if these magnetic fields have a large extent (L 3> Mpc) and a significant 
strength 10~^ G), their effects could produce noticeable signals at energies above 
10^° eV even for protons, while for more modest fields, it would be the galactic 
magnetic fields the ones producing the dominant effects. 
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